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Abstract 

We show that the standard perfect fluid paradigm is not necessarily a valid description of a 
steady state gravitational source unless the background geometry has very high symmetry. By 
incoherently averaging over a compete set of modes of a scalar field propagating in a curved 
background, we show that for a static, spherically symmetric geometry the energy-momentum 
tensor that ensues will in general be of the form T^y = {p + p)U^Uu + pg^u + QVfiVu where 
is a spacelike vector. Such a qV^Vy type term is absent for an incoherent fiuid in flat spacetime 
since a Minkowski geometry has much higher symmetry than a spherically symmetric one, and 
would thus be missed in a covariantizing of a fiat spacetime T^jy. However, in a general relativistic 
treatment of systems such as stars, such qV^Vu type terms cannot automatically be ignored, with 
their contribution in the strong gravity limit potentially being quite significant. 
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I. INTRODUCTION 



In gravitational theory it is standard to take the energy-momentum tensor of a steady 
state gravitational fluid source to be in the form of a perfect fluid. As such, a steady state 
fluid consists of a set of particles which are all in free fall or a set of wave modes whose 
contributions to the energy-momentum tensor are added incoherently. In flat spacetime the 
energy-momentum tensor of such a fluid will have the perfect fluid form 

Tf,u = {pf + Pf)Uf,U^ + PfQf^u, (1) 

with fluid energy density pf and fluid pressure pf, and it is conventional to covariantize this 
expression in order to obtain the energy-momentum tensor which is to describe a steady 
state fluid in a curved geometry. The covariantization prescription is to replace all flat 
space tensors by curved space ones and all flat space derivatives by covariant ones, so that 
a curved space perfect fluid is to be described by ([1]) as written in a non-flat background, 
with no other terms with the requisite tensor structure being deemed relevant. Remarkably, 
no justification for the appropriateness of such a procedure appears to have been given in 
the literature, and its use is simply taken as being self-evident. 

In this paper we examine this prescription and find that it does not in fact hold in 
general. In particular, we find that it does not necessarily hold in a situation where it is 
commonly used, namely in a spherically symmetric system such as a star. To understand 
why there might even be a concern for spherical systems, we note that spherical systems 
are only isotropic about a single point (the center of the system). Such a sjTumetry only 
requires for any given tensor A^^, that its Ag^ and A^'^ components be equal, and imposes 
no restriction on A/. A familiar example of this is the form of the Einstein tensor G^,^ in 
the static, spherical geometry 

ds^ = -B{r)df + A{r)dr^ + r^dO^ + sin^ Odcj)^, (2) 

where the non-zero components of G^^p are given by 
^ _ B B BA' 

Grr = - 

_ G^^ _ r^B" r^A'B' r^B'^ rB' rA' 

~ ih?^ ~ ~2AB ^ T¥b ^ Ta& ~ 2AB ^ 2A^' ^ ^ 
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As we see, there no relation of the form G/ = Gg^ in Q, and not only that, there could 
not be since the Bianchi identity G'^'^.^ = relates the radial derivative of G"''" to G^^, to 
thereby force G^^ to be a first-derivative function of the metric coefficients and to be a 
second-derivative one. This same problem is not evaded if one works in isotropic coordinates, 
where, because of coordinate invariance, the metric can actually be brought to a form 

ds'^ = -H{p)dt^ + J(p) (dp^ + p'^de'^ + sin^ Odcf)'^) (4) 

that does have the generic perfect fluid form g^P = g^^ = g^'^. Specifically, even in this 
coordinate system the Einstein tensor 

2HJ' 3HJ'^ HJ" 



Gqo — — 77; r~r^ — I" 

G 



pp 4J3 J2 ' 

- H'J' J'^ H' 

'''' ~ ^ 2EJ ~^~~pE 



2J 2J 2J2 2H Am 2H ^ ' 

is still not of a perfect fluid form, and indeed must still not be since now the Bianchi identity 
relates the radial derivative of G'''' to G^^ . Since the gravitational equations of motion 
would relate the Einstein tensor (or some equally covariant alternate gravity equivalent) 
to the energy-momentum tensor, even in isotropic coordinates we should not expect the 
components of T^y to necessarily obey T^^ = Tg^ = T^'^ . 

To determine what should be the case, and to understand why one is tempted to use ([1]) 
in the curved space case, it is instructive to ask how ([1]) came about in the flat space case in 
the flrst place. Speciflcally, to establish a flat space perfect fluid form, one has to evaluate 
the partition function for the system, and to do so one incoherently averages over all of the 
modes of the system. When this is done for the modes of a scalar field or of a Maxwell field 



propagating in a flat background, one does indeed recover the perfect fluid form [l, 2], just 
as one should since it is precisely such an incoherent sum over modes which produces the 
energy and pressure of flat space black-body radiation and the familiar pf = Pf/3 associated 
with a traceless T^jy. The interest of Q and Q (and also [3!]) was to extend the incoherent 
averaging calculation to a background Robertson- Walker geometry, with ([T]) indeed then 
being recovered, just as one would want. However, both a Minkowski background and a 
Robertson- Walker background have very high symmetry, Robert son- Walker being maximally 
3-symmetric and Minkowski being maximally 4-symmetric, with both of these geometries 



being isotropic about every spatial point. In contrast, unlike these particular geometries, 
the geometry associated with ([2]) or (jl]) has much lower symmetry, being isotropic only 
about one point, the origin of coordinates. We thus cannot immediately anticipate that an 
incoherent averaging over modes propagating in the geometry associated with ([2]) or (jlj) will 
necessarily lead to ([T]). In fact, we will show that in general it does not necessarily do so. 
Specifically, to show that we do not always recover ([1]) in general, we only need to find one 
counterexample. We thus only need to find one appropriate choice of metric coefficients for 
which one can do the incoherent summation over an infinite set of modes analytically and 
then fail to recover ([T]). It is precisely such a calculation that we provide in this paper. 

In Sec. HI] we introduce the incoherent averaging procedure. In Sec. IIIII we apply it to a 
fluid in a specific, exactly solvable, static, spherically symmetric geometry, to find that we 
do not obtain a perfect fluid form. In Sec. IIVI we discuss the difference between treating 
fluids as incoherently averaged waves or treating fluids as incoherently averaged particles 
that are in geodesic free fall. In Sec. |V] we show that through the imposition of boundary 
conditions there can be departures from a perfect fluid form for finite-sized systems even 
in flat spacetime. Finally, in Sec. |Vl] we discuss the implications of kinetic theory for the 
structure of fluids. 



II. THE CURVED SPACE ENERGY-MOMENTUM TENSOR 

For our purposes here an appropriate model to study is a minimally coupled scalar field 
S{x) with action 

I = -Jd'xi-gy/'^V,SWS (6) 

propagating in the background associated with the metric of (jlj). For the scalar field the 
equation of motion is given by V^V^S = 0, i.e. by 

_ l_d^ 1 d 

(7) 

and the energy-momentum tensor is of the form 

T,. = V,SV,S - ^(7m.V,^V"5. (8) 



op 



+ 



Jp^ 



+ 



1 d^S 
sin^ 9 dctP' 
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The equation of motion can be separated, and for a real scalar field the general solution can 
be expressed in terms of four characteristic solutions 

Si{x) = NP sm{ujt)S^^g{p)PP{e) sin(m0), ^2(2;) = sm{ujt)S^^e{p)PP{e) cos(m0), 
Ss{x) = Np cos{ujt)S^^e{p)Pl^{e) sin(m0), S^{x) = cos{ujt)S^^e{p)Pp{e) cos(m0), 

(9) 

where N]!^ is given by 



{2i+l){i-m)\ 



and where the radial term obeys 



1/2 



+ 1) 



Jp2 



Su;,i{p) = 0. 



(10) 



(11) 



To illustrate the incoherent averaging procedure, we note that in a given mode such as 
5*1(0;) a quantity such as VaSV'^S evaluates to 



^2 1 
H J 



dSu. 



1 2 



dp 



[NpppYsm\m<, 



+ 



1 



sin^(a;t)[S'^,^ 



mi 2 



dPl 

de 



m 



sin^(m0) + — ^[Nppp]\os\m(, 



.(12) 



On now adding to this expression those obtained when one uses the other three above 
solutions, one obtains 



i=4 



i=i ^ 



dS„ 



Jp- 



mi 2 



dPl 

de 



dp 
2 



m omi2 



[NTP, 



m 



m Dml2 



sin'^ 



(13) 



Using standard properties of the spherical harmonics, 

(2£ + 1) 



mi2 r D"^l 2 



47r ' 
(2£ + l)£(£ + l)sin2^ 



57r 



ml 2 



dPl 

de 



(2£ + l)£(£ + l) 



57r 



(14) 



just as in [11] one sums over the azimuthal quantum number m to obtain 

2" 



m i=l 



An 



1 
J 



dS, 



dp 



Kiu,p,i), (15) 
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with (fTS!) serving to define K{lj, p, £). With regard to (fTSi) . we note that because the angular 
part of the metric is maximally 2-symmetric, the sum on the azimuthal m removes any 



dependence on the ang^ 
yields (again following 

Tpp{uj,pJ) 

Tgg{uJ,p,i) 



e 9. Repeating this same procedure for the rest of T^^, of 

HK 



then 



(2£ + 1) 



+ 



dS, 



1 2 



dp 



2 ' 
JK 



sin' 9 



T^^{uj,p,i) 



2 ' 
(2£+l 



+ 1) 



Stt 



p\JK 



(16) 



with all other components of T^^, being zero. 

To complete the incoherent averaging we need to sum over all i values as well, an infinite 
summation. However, if we revert back to flat space where H{p) = J{p) = 1, the radial 
equation is then solved by the spherical Bessel functions ji{ujp), and the sum on £ can then 
be performed analytically using the completeness relations for Bessel functions: 



f dp 



' dp^ 



UJ 



E(2^ + i) 

2cuV' 



dp 



UJ 

y 



(17) 



On defining Kiyj^ p) = J2e K{uj, p, i) and T^u{uj, p) = J2e T^uiu, p, i), one obtains K{uj, p) = 
and 



Too{uj,p) 



Tpp{uJ,p) 



1 uo'^p^ 
Teg{u}, p) = —YT^Tni^, p) = -^T^, (18) 



47r' i2vr' -^v-"- sin^ ' 127r 

One thus obtains none other than the perfect fluid form given in ([1]), with the averaging on £ 
removing any dependence on the coordinate p from T^" because of the maximal 3-symmetry 
of the spatial part of a Minkowski metric. The incoherent averaging prescription thus gives 
a perfect fluid in flat spacetime, just as one would want. 



III. INCOHERENT AVERAGING IN CURVED SPACETIME 

To follow this same procedure in curved space is not at all as straightforward, since for 
an arbitrary choice of H{p) and J(p) the radial equation will not necessarily be solvable in 
terms of named functions, and the needed completeness relation for the ensuing modes may 
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not even be known at all. Moreover, in the dynamical case where the energy-momentum 
tensor is used as the source of tensors such as the Einstein tensor given in ([5]), one has 
to solve for H{p) and J(p) self-consistently, a procedure that would not only not appear 
at all likely to yield an analytic result, since it involves an infinite summation, it does not 
immediately lend itself to numerical approximation either. However, to test the viability 
of the perfect fluid assumption, one only has to seek an appropriate choice of H{p) and 
J{p) for which one can do the £ summation analytically, to see whether or not the relation 
Tpp = Tqq I then does in fact ensue. 

Since we are able to do the I summation analytically for Bessel functions, we shall seek a 
form for -ff (p) and J(p) for which the radial equation can be reduced to the Bessel equation. 
To this end we set J(p) = Hi^p) = P{p)- For this choice, the metric in (jl]) reduces to 

ds" = f{p)[-dt^ + dp" + p'de^ + p2 sin' ed(j)% (19) 

to thus be conformal to flat. However, it is not flat since in it the Einstein tensor in ([5]) 
takes the non- vanishing form 



00 



4/' f" ^ 2f" 



G 



4f 3/'^ 

Pf 



Gf) 



2pf ^ p'P 



2pV" 
/ 



pf P ' / ' pf f ' / ' P 

a form which is still not a perfect fluid unless /(p) just happens to obey 

/ pf P 



0. 



(20) 



(21) 



For the metric of (fT9l) . on substituting S^j{p) = Qui/ip)/ f{p) the radial equation (fTTj) 
reduces to 



d^ 



d 



~^ dp^ ~^ dp 
Thus, for any choice of /(p) for which 



+ 1) 1 



2/' 



P' 



-Af" + — 
f \ p 



QuAp) = 0- 



(22) 



(23) 



where is positive, one finds that (!22|) reduces to none other than the Bessel function 
equation 

r _ _ ^2 9 ^ ()(() _i_ 1 \- 



cj' + k' + 



dp"^ ^ pdp 



p- 



with immediate solution Quj,i{p) = jii^p) where A 
readily given as 

fip) 



QluAp) = 0) 

uJ^ + K^y/\ For 



[asin(fi;p) + /3cos(Kp)] 



P 



solutions are 



(25) 



and thus form a whole family of solutions labeled by all real values of the parameters k, a 
and /5. For any /(p) which obeys fl23|) . the Einstein tensor reduces to 
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f 

Goo = — — — 2k 



__4f _3r 
Pf P - 



Goo = ^ + ^^ + 2p\', (26) 



P ' pf P ' f f 

and in solutions of the form given in ( l25l) is still not in the form of a perfect fluid, with the 
solutions in fl25|) not obeying fl2T|) . 

To now evaluate the energy-momentum tensor when (123|) is imposed, on summing over £ 
as before, one obtains 

K{u^,P) = K^^ (27) 



4vr/4 ' 

with the incoherently-averaged energy-momentum tensor itself being given by 



1 


+ + 


k2 ■ 


47r 


2/2 


1 


' 

3/2 ^ 2/4 


k2 


Air 


6/2 



^ee(t^,p) 



1 

47r 



2 fl2 



P'f 



^2/^2 



(28) 



3/2 2/4 6/2 

As we see, T^^{uj,p) is not in the form of a perfect fluid since T^f — Tg^ = f^/Anf^ is 
not zero. Since our analysis is valid for any /(p) which obeys ( l25i) . we recognize a whole 
family of metrics for which a perfect fluid is not obtained. In general then, the use of perfect 
fluid sources for spherically symmetric gravitational systems must be regarded as open to 
question. 



IV. CONTRASTING AVERAGING OVER MODES WITH AVERAGING OVER 
PARTICLES 

To characterize the tensor structure of the T^^, that incoherent averaging has led us to, 
we need to augment the perfect fluid energy-momentum tensor T^j, of ([1]) with a further 
term of the form qjV^Vy where is a unit spacelike vector in the radial direction, and just 
like Pf and p/, qj is equally a general coordinate scalar. With reference to the metric of (jll), 
we define = (0, and with f/^ = {H^^"^, 0,0,0) as usual, find that we can write 

the energy-momentum tensor as 

Tf,, = {pf + Pf)U^U, + pfQ^, + g/ V^K- (29) 
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For the energy-momentum tensor of (l28l) we identify 



1 

Pf 



Arc 



2/6 2/ 



1 



47r 



3/4 2/6 6/4 



f/2 



with a summation over uj being imphcit. For the form of T^i, given in (129|) . in the geometry 

associated with the metric of (jlj) the covariant conservation condition T^'^.^ = yields 

(jf 2 \ ff' 
7 + - j + ^(P/ + P/ + (If) = 0, (31) 

to thus show the exphcit role played by the qf term in the energy and momentum balance. 

Not only does qj act analogously to pf in the covariant conservation condition, it even 
does so in dynamical equations such as the Einstein equations G^^, = —SnGT^^. Specifically, 
if we work in the weak gravity limit where we can set H{p) = 1 + h{p), J{p) = 1 + j{p) 
with h{p) and j(p) both being of order G, we can consistently find solutions in which pf 
is of order one and pf and g/ are both of order G. In the weak gravity limit the Einstein 
equations and the conservation condition associated with the metric of (jlj) reduce to 

— + / = -SttGp/, - - = -MG{pf + qf), 

P P P 

?■" ?■' h' h" 2 h' 

-J-Tp-Tp-Y- + + -p^' + 2 = °' 

and to order G thus have a consistent weak gravity limit in which j + h = 0. In and of 
themselves then, the dynamical equations do not require g/ to be at least one order in G 
smaller than pf. 

In discussing fiuids one can consider them to be composed of either particles or waves. To 
effect the incoherent averaging in the above we averaged over waves rather than particles. 
However, for a set of non-interacting particles, if one takes them to be in free fall and 
averages over them, one does get a perfect fiuid with no qf type term. Specifically, since 
no such term is present for such a fiuid in fiat spacetime where all the particles move on 
fiat space geodesies, in the presence of gravity, free-falling particles will continue to move 
on geodesies and thus still behave as a perfect fiuid. We thus need to ask why it was that 
the averaging over waves did not give a perfect fiuid, and how one can then reconcile the 
free-falling particle and the incoherent-averaged wave descriptions. To this end we need to 
recall the conditions under which rays normal to wavefronts are geodesic. 

n 

As noted in |4[], if we set S{x) = exp{iT{x)), from the wave equation V^V^S = we 
obtain 

VTV^T - iVV^T = 0. (33) 



For wavelengths that are short with respect to the typical scale of the problem, (l33l) reduces 
to V^TV^T = 0, to yield first V^TV^V^T = 0, and then 

V^TV^V^T = (34) 

since V^V^T = V^'V^T. Since normals to the wavefronts obey the eikonal relation 

WT = ^ = F (35) 
dq 

where q measures distance along the normal and /c^ is the wave vector of the wave, on noting 
that d/dq = {dx^/dq){d/dx^), from fl34l) we thus obtain 

^"V„.' - ^ + - 0. (36) 

Recognizing (l36l) as the massless geodesic equation, we see that in the short wavelength 
limit, rays move on geodesies. Moreover, exactly the same result can be obtained jj] for 
massive scalar fields as well, with it being a general rule that short wavelength rays are 
geodesic. Referring now to the incoherently averaged T^^, given in fl28l) . we see that in the 
large u limit T^,y does indeed reduce to a perfect fiuid with the g/ term becoming negligible. 
Perfect fiuids are thus to be associated with short wavelength modes alone, with departures 
from a perfect fiuid form being expected to occur at longer wavelengths where one can no 
longer ignore the V^V^T term in ( l33l) . The potential importance of such long wavelength 
modes thus depends on the scale of spatial variation of the system of interest, a dynamical 
rather than a kinematical issue. Moreover, the relative importance of the longer wavelength 
modes will increase as the size of a system is decreased, to potentially give a departure from 
a perfect fiuid form that could be substantial when a system such as star collapses to a size 
of order its Schwarzschild radius j^. 



V. IMPLICATIONS OF BOUNDARY CONDITIONS 

While we have discussed possible departures from a perfect fiuid form in curved space, 
such departures can even occur in fiat space due to boundary effects of finite sized systems. 
Such concerns do not arise in fiat space when we consider plane waves, since they fill all space. 
However, suppose we consider a finite-sized cavity which is filled with massless scalar modes 
in thermal equilibrium (a spinless black body). Now it is long known that if we take the 
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cavity to be a cubical cavity of side of length a and take the modes to obey periodic boundary 
conditions, while we would find corrections of order hc/kTa to the familiar Too ~ T'^ black- 
body formula, because of the cubic symmetry these corrections would treat all components 
of the pressure tensor equivalently and leave the perfect fluid form T^^ = Tyy = T^z = Too/3 
(i.e. Trr = Tog/r"^ = T^fjj/r'^ iAv? 6 = Too/3) intact. These finite size effects would only be 
significant for wavelengths of order a, and would become inconsequential as the size of the 
cavity is increased. 

However, suppose we consider a spherical cavity of radius a, and rather than periodic 
boundary conditions Gj, instead require that the modes have radial wave functions which 
vanish at the surface of the cavity. The radial functions would then obey je{uja) = 0, with 
the allowed frequencies then being given by the zeroes of the Bessel functions. Each Bessel 
function has its own infinite set of discrete zeroes [e.g. the zeroes of jo{x) = (sina;)/x obey 
sinx = 0, and those of ji{x) = (sinx)/x^ — (cosx)/x obey sinx = xcosx]. On labeling 
these zeroes as j^, inside a flat space spherical cavity the allowed frequencies are given by 
the discrete a;.^ = ji/a, and the radial wave functions themselves are given by je{jir/a). [In 
fiat space there is no distinction between the radial coordinates r and p of ([2]) and (jlj).] The 
sum over modes proceeds initially as previously, with f|T6|) being replaced by 



(2£+l) 
Stt 

(2£ + 1) 
Stt 

(2^ + 1) 
Stt 



J\2 



+ 



( l\2 



+ 1)' 



/ c(7^(jn^/Q) y 
I dr ] 



dr 

' djeijjr/ay 
dr 



(37) 



To complete the incoherent averaging we would need to sum over all uj^ for a fixed i and 
then sum over all i. However, to see whether or not a perfect fluid form might emerge, we 
only need to look near r = a. Noting that every Bessel function obeys the recursion relation 

djeiujr) 



dr 



r 



ujje+i[ujr), 



(38) 



and recalling that the zeroes of jilx) are distinct from those of j£^i{x), we see that at r = a 
jiiJir/a) vanishes but dji{j^r /a)/dr does not. Thus at r = a we have 



(2£+l) fdUjir/a) 



in 



dr 



(39) 
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FIG. 1: Plot of T„ (dashed curve) and TqqIt^ (continuous curve) as a function of r as incoherently 
summed over the hrst 100 zeroes of each of the first 100 I values of jeiju^/a) in a spherical cavity 
of radius a = 1. 

Since Trr{uj^,a,i) is positive definite and equal and opposite to the negative definite 
T0g{ul^,a,£)/a^, and since the sum over c<j^ and £ does not change this, we see that we 
do not recover the perfect fluid form at r = a. A straightforward Taylor series expansion of 
the form r = a — e shows that this result remains true to second order in e. To get a sense 
of how the incoherent sum at arbitrary r might look, we have numerically summed over the 
flrst 100 zeroes of each of the flrst 100 i values of ji{jir/a) (i.e. 10,000 terms in total), and 
as we see in Fig. (1), there is no sign of a perfect fluid form. Even in flat spacetime then, 
boundary conditions can prevent a fluid from being of the prefect fluid form when the fluid 
is in thermal equilibrium inside a spherical cavity of flnite radius. 
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VI. IMPLICATIONS OF KINETIC THEORY 



While the gravitational equations themselves provide no specific basis for leaving out 
any type term from the fluid T^^ in general, in the weak gravity limit one might still 
be able to exclude such terms via the use of kinetic theory. In applying kinetic theory to 
gravitational systems there are two approaches, one based on the Boltzmann equation, and 

the other based on the Liouville equation. Since there are some differences between these 

i 

two approaches [?]] with the former being more appropriate for stars and the latter more 
appropriate for systems such as a cluster of galaxies, we describe them both. 

Consider first a set of particles each of mass m in some long-range (typically gravita- 
tional) external potential Vext{^) that are undergoing rapid (typically atomic) momentum 
conserving two-body collisions v + V2 ^ v' + V2 through a scattering angle Q with differ- 
ential collision cross section a{Q). In the absence of two-body correlations the one-particle 
distribution function /(x, v,t) obeys the Boltzmann equation (see e.g. 8j) 

g/(x,v,t) ^ ^ g/(x,v,t) gKxf(x) g/(x,v,t) 

dt (9x 9x (9v 

= / |v-V2|a(r])[/(x,v,'t)/(x,v^,t)-/(x,v,t)/(x,V2,t)]rf3v2rfa (40) 

As such, fl40p admits of an exact time-independent Maxwell-Boltzmann type solution 



/mb(x, p,)!:) = Cexp 



mv^ mVn^t fx) 



(41) 



2kT kT 

where the temperature T and the coefficient C are independent of x. (That /mb is an exact 
solution is because it makes both sides of (HOl) vanish separately.) However, since averages 
with this distribution function are given by (A) = J d^pAfuB/ J c^^p/mb, for any dynamical 
variable y4(x, v,t), the Kxt(x) term would drop out and all averages evaluated with this 
/mb would be spatially independent. Moreover, the number density itself would behave 
as r2(x) = / (i^p/MB(x, p, t) ~ exp(— mVext(x)/A;T), and give a spatial dependence that is 
nothing like that required of the density of a star of finite size. 

Now as such, use of the above Boltzmann equation presupposes that the only collisions 
of relevance are atomic type ones, and that those collisions dominate the relaxation of the 
system to the distribution function /mb(x, p, t). However, in stars gravitational collisions 
also play a role, and their effect cannot be isolated solely in the c^Kxt (x) /(9x term in (HOll . 
Rather, they serve to modify the right-hand side of (HOj) as well. However, since the cross- 
section for gravitational scattering is infinite, the effect of gravity cannot be accounted for by 
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a collision integral term in which one simply includes a gravitational contribution to (t(Q). 
Rather, one should work with the Liouville equation as it is better suited to handle long 
range forces, and we will discuss this below. 

However, before doing so, we note that in kinetic theory, through use of the method of 
the most probable distribution, it is possible to determine the distribution function without 
needing to actually solve or even construct the Boltzmann equation at all. This method does 
not determine the approach to equilibrium (i.e. the temporal behavior of the distribution 
function), but does give the equilibrium configuration that results, and it is valid even if the 
distribution function does not obey an equation such as ( HOj) at all. Thus we can use the 
method of the most probable distribution to get a sense of how gravitational interactions 
might modify the distribution function fHT]) . Ordinarily one applies the method by taking 
the system of interest to be confined to a region of phase space with a fixed total energy 
and fixed volume and introduces spatially-independent Lagrange multipliers for the total 
number of particles and the total energy j^. In the absence of gravity this leads to the 
Maxwell-Boltzmann distribution being the overwhelmingly likely one. In the presence of 
non-relativistic gravity we can adapt this method to concentric shells within a star, and 
take the temperature and density to be a constant within any given shell but to vary from 
one shell to the next. The approximation here is that particles stay within shells and do not 
exchange energy or momentum with any particles except those in their own shell. (Since 
gravity is a long range force, at best this assumption could only be valid when gravity is 
weak.) In this approximation the Lagrange multipliers for the total number of particles 
and the total energy within a given shell are taken to be uniform throughout the shell but 
to depend on the location of the shell, with the most probable (mp) distribution function 
throughout the star then being given by 

\2" 



/""P(x, p, t) = — —j^ exp 



m V — u 



29 



(42) 



In f H2|) the particle number density n(x, t) = / (i^p/™P, the average velocity u(x, t) = 
/ (i^p/™Pp/mn(x, t) and the temperature 6'(x, t) = (m/3) / (i^p/™P|v — ^(x, t)p/n(x, t) are 
now all spatially dependent. 

While such a distribution function would not be an exact solution to the Boltzmann equa- 
tion, in the event that n(x, t), u(x, t) and 6'(x, t) are all slowly spatially varying, /™p(x, p, t) 
would be a good first-order approximation to it. For such a distribution function the pressure 
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tensor evaluates to 



Pij = mn{x, t){{vi - Ui){vj - Uj)) = 5ijn{x, t)9{x, t), (43) 

to thus be of none other than of perfect fluid form. If the full distribution function is taken 
to obey the Boltzmann equation, then in the next order the pressure tensor evaluates to 81] 



Pij = 5ijn9 — rnO 



dui ^ duj 2^ ^ 
dxj dxi 3 



(44) 



where r is the mean free time between particle collisions. With r;,(x, t) depending on position, 
the djUi + diUj term would not be proportional to 5ij. When n(x, t) is slowly varying then, 
the pressure only begins to depart from 5ij in corrections to first order. To the extent then 
that weak gravity would cause a star to bind with a slowly varying density and that r would 
be small (i.e. a small mean free time between atomic collisions), the g/ type term could be 
neglected in lowest order. Thus our ability to ignore g/ in a weak gravity star depends on 
how good a dynamical approximation the above /™p(x, p, t) is to the full /(x, p, t). The use 
of perfect fluids as gravitational sources for weak gravity stars is thus equivalent to using 
/™P(x, p, t) as the one-particle distribution function. However, once one has to go beyond 
/™P(x, p, t) as would be the case for stronger gravity, there would no longer appear to be 
any immediately apparent justification for ignoring qj type terms. 

In terms of the geodesic equation discussion given earlier, we now see that short wave- 
length for eikonal purposes means short with respect to the distance scale on which the 
particle number density n(x, t) varies. The slower the spatial variation of n(x, t) then, the 
fewer the number of modes that will not incoherently average to a perfect fluid. 

In the above discussion it is the non-gravitational collisions which dominate, with the 
discussion being given for a weak gravity star in which there are atomic collisions between the 
atoms in the star. However, for a cluster of galaxies, it is gravity itself which has to establish 
an equilibrium distribution of galaxies. In clusters there is a lot of X-ray producing plasma 
located in the region between the individual galaxies in the cluster. Collisions between the 
atomic particles in the plasma can readily bring the plasma to thermal equilibrium, but 
unless they can bring the galaxy distribution to equilibrium too, it would be up to gravity 
to do so. 

To describe the role that gravity plays in a cluster of galaxies we turn to the Liouville 
equation. We treat each of the N galaxies in the cluster as a non-relativistic point particle 
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of mass m with position Xj and velocity Vj. Each galaxy moves in the gravitational field 
0(xj) produced by the other galaxies in the cluster, and obeys 

One introduces the normalized (to one) 6A^-dimensional distribution function 
/'•^•'(xi, vi, ,XAr, VAr,t), and finds that it obeys the Liouville equation 



0, 



(46) 



where 

(Pai^a) = 0(Xa,X^). (47) 

If the distribution function is symmetric under the interchange of any pair of particles and 
sufficiently convergent asymptotically, upon integrating (146!) . one finds (see e.g. [Qj) that for 
large N, the one- and two-particle distribution functions 



/«(xi,Vi,t) = J d'x2dW-ANd\Nf^''\ 

/(2)(xi,Vi,X2,V2,t) = J d'x,dW..ANd^^^Nf^''\ 



are related by 



+ vi ■ — — = N 



dt 9xi J dxi dvi 

In terms of the two-body correlation function defined by 



(48) 



(49) 



5((Xi, Vi,X2, V2,it) = /( )(Xi, Vi,X2, V2,t) - f ^ (Xi , Vi , t) )(X2, V2,)f:) 

and the kinetic theory distribution /(xi, vi,t) = Nf^^\xi,\'i,t) that is normalized to 

J (i^v/(x, V, t) = n{x, t), J d^xn{x, t) = N, 



(50) 



we find that fj49l) takes the form 

9/(x,v,t) 



dt 



+ V- 



9/(x,v,t) dV{x) 9/(x,v,t) 



9x 



dv 



av 



d f d(j){x,X2] 



(9x 



5((X, V,X2, V2,t)rf^X20?V2 



where 



V{x) = J C?^X2rf^V2/(x2, V2,t)0(x,X2) = J rf^X2n(x2 , t)0(x, X2) . 



(51) 



(52) 



(53) 
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As such, the potential V^(x) introduced in (l53l) serves as a mean-field potential, and should 
the system relax so that the two-body correlation function becomes negligible, the left-hand 
side of ( 1321) would then become equal to zero. At first glance the left-hand side of ( 13^ 
looks quite like the left-hand side of the Boltzmann equation fHUl) . However, in fHUl) \4xt(x) 
is an external one-body potential, while in (132]) V(x.) is a statistically averaged two-body 
potential. Moreover, while the vanishing of 5'(xi, vi, X2, V2, t) would cause the right-hand side 
of fl52|) to vanish, it would not have the same effect on the right-hand side of the Boltzmann 
equation of (l40l) . In fact, it was the very requirement that 5'(xi, vi, X2, V2, t) vanish that led 
81] to the explicit form for the collision integral term given on the right-hand side of ( HOj) 
in the first place, with the vanishing of 5'(xi, Vi, X2, V2, t) not requiring the vanishing of the 
Boltzmann equation collision integral term. When the right-hand side of fl32]) does vanish, 
all one can say is that the steady state solution to fl52l) simply requires that /(x, v,t) be 
an arbitrary function of the quantity -|- V(x). With (1521) possessing no analog of the 
collision integral term in f j40l) (gravity being l ong range), there is nothing to force /(x, v,t) 
to be an exponential function of v^/2 -|- V{x.) [10]. The specific dependence on + V^(x) 
that /(x, v,t) would acquire in equilibrium would depend entirely on how the two-body 
(^(xi, Vi, X2, V2, t) would behave before it becomes negligible. In the absence of knowing 
how 5'(xi, vi, X2, V2, t) behaves, one is not able to recover any analog of (H3l) . Consequently, 
without detailed information on the two-body correlation function, for clusters, even weak 
gravity ones, one is not able to anticipate that in steady state the galaxies in a cluster would 
behave like a perfect fiuid 11]. 
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